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Abstract 

In this paper, we study a multi-dimensional backward stochastic differential equa- 
tion (BSDE) with oblique reflection, which is a BSDE reflected on the boundary of a 
special unbounded convex domain along an oblique direction, and which arises natu- 
rally in the study of optimal switching problem. The existence of the adapted solution 
is obtained by the penalization method, the monotone convergence, and the a priori 
estimations. The uniqueness is obtained by a verification method (the first component 
of any adapted solution is shown to be the vector value of a switching problem for 
BSDEs). As applications, we apply the above results to solve the optimal switching 
problem for stochastic differential equations of functional type, and we give also a prob- 
abilistic interpretation of the viscosity solution to a system of variational inequalities. 
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1 Introduction 



In this paper, we are concerned with the following reflected backward stochastic dif- 
ferential equation (RBSDE for short) with oblique reflection: for i 6 A := {1, • • • , m}, 
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Yi(t) = ii + f ^(s,Y i (s),Z i (s),i)ds- J dKi(s)- J Zi(s)dW(s), t€ [0,T], 
< Yi(t) < mm{Y;-(t) + fc(*,j)}, 

[ (y^s) - mm{Yj{s) + k(i,j)}) dK z (s) = 0. 

(1.1) 

Here, W is a standard Brownian motion on a complete probability space (Q, J 7 , P) and 
£ is an m-dimensional random variable measurable with respect to the past of W up 
to time T. £ is the terminal condition and ip the coefficient (also called the generator). 
k is a real function defined on A x A. The unknowns are the processes {^(£)}te[o,T] , 
{Z(t)} te i 0jT i, and {K(t)} te i 0T i, which are required to be adapted with respect to the 
natural completed filtration of the Brownian motion W. Moreover, K is an increasing 
process. The third relation in (1.1) is called the minimal boundary condition. 

RBSDE (1.1) evolves in the closure Q of domain Q: 

Q ■= {(yi, ■■■ , Vmf eR m :yi< yj + k(i,j) for any i, j G A such that j ^ i}, 

which is convex and unbounded. The boundary dQ of domain Q consists of the bound- 
aries dL~,i G A, with 

L 7 '■= {(yi> ■ • • > Vmf eR m :yi< yj + k(i,j), for any j G A such that j A. 
That is, 

m 

dQ = U 0L" 

In the interior of Q, each equation in (1.1) is independent of others. On the boundary, 
say dL~, the i-th equation is switched to another one, and the solution is reflected 
along the oblique direction — a (which is the opposite direction of the i-th coordinate 
axis) . 

RBSDE was first studied by El Karoui et al. [9] for the one-dimensional case. Mul- 
tidimensional RBSDE was studied by Gegout-Petit and Pardoux [10], but their BSDE 
is reflected on the boundary of a convex domain along the inward normal direction, 
and their method depends heavily on the properties of this inward normal reflection 
(see (l)-(3) in [10]). We note that in a very special case (e.g., tp is independent of z), 
Ramasubramanian [16] studied a BSDE in an orthant with oblique reflection. Note 
also that there are some papers dealing with SDEs with oblique reflection (see, e.g. 
[13] and [7]). 

An incomplete and less general form of RBSDE (1.1) (where the minimal condition 
of (1.1) is missing and the generator ip does not depend on (y,z)) is suggested by [3]. 
But they did not discuss the existence and uniqueness of solution, which is considered 
to be difficult. See Remark 3.1 in [3]. 

Besides the theoretic interest, RBSDE (1.1) arises naturally from the following 
optimal switching problem. 
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Consider the switched equation 

X<\t) = x + f a(s,X a ^)[dW(s) + b(s,X<\a(s))ds], t e [0,T] (1.2) 
Jo 

and the cost functional 

J(a(-))=E[[ l(s,X a (-\a(s))ds]+E[Y / k(a i - 1 ,a i )}. (1.3) 

The optimal switching problem is to minimize the cost J(a(-)) with respect to a(-), 
subject to the state equation (1.2). 

In the above, xo is a fixed point in R d . a, b and I are defined on [0, T] x C([0, T};R d ), 
[0,T] x C([0,T] ;j R d ) x A and [0, T] x C([0,T];i? d ) x A, respectively, with appropriate 
dimensions. 

oo 

a(0 = "oX{e„}(-) + 5^ai-iX(e i _i,9 i ](0 
i=i 

is called an admissible switching strategy if for any i, 9i is a stopping time, and on is 
an J T g i -measurable random variable with values in A. Here, x is the indicator function. 
k is called the switching cost. 

Optimal switching is a special case of impulse control. The classical method of 
quasi-variational inequalities to solve impulse control problems driven by Markov pro- 
cesses is referred to the book of Bensoussan and Lions [2]. See [18] and the references 
therein for the theory of variational inequalities and the dynamic programming for op- 
timal stochastic switching. But these works are restricted within the Markovian case. 
Recently, using the method of Snell envelope (see, e.g. El Karoui [8]) combined with 
the theory of scalar valued RBSDEs, Hamadene and Jeanblanc [11] studied the switch- 
ing problem with two modes (i.e., m = 2) in the non-Markovian context. Djehiche, 
Hamadene and Popier [6] generalized their result to the above switching problem with 
multi modes. We note that in both [11] and [6], the drift term b does not depend on 
the switching strategy a. 

The main contribution of this paper is to establish the existence and uniqueness 
of solution for RBSDE (1.1). We prove the existence by the penalization method, the 
monotone convergence, and the a priori estimation whose proof is rather technical. The 
proof of uniqueness is quite different: the classical method to estimate the difference 
of two solutions appears difficult to be applied to our present case of the oblique 
reflection. We obtain the uniqueness by a verification method: first we introduce an 
optimal switching problem for BSDEs, then we prove that the first component Y of 
any adapted solution (Y, Z, K) of RBSDE (1.1) is the (vector) value for the switching 
problem. As applications, we solve the optimal switching problem (1.2) and (1.3), 
and we establish the Feynman-Kac formula for the viscosity solution to a system of 
variational inequalities. 

The paper is organized as follows: in Section 2, we prove the existence of solution, 
whereas Section 3 is devoted to the uniqueness. We solve the optimal switching problem 
in Section 4. Finally, in Section 5, we give a probabilistic interpretation of the viscosity 
solution to a system of variational inequalities. 
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2 Existence 



2.1 Notations 

Let us fix a nonnegative real number T > 0. First of all, W = {Wt}t>o is a stan- 
dard Brownian motion with values in R d denned on some complete probability space 
(£l,J-,P). {J~t,t > 0} is the natural filtration of the Brownian motion W augmented 
by the P-null sets of T . All the measurability notions will refer to this filtration. In 
particular, the sigma-field of predictable subsets of [0, T] x Q, is denoted by V. 

Let us consider now the RBSDE (1.1). The generator tp is a random function 
■ip : [0,T] xflxfix R d x A — > R whose component ip(-, i) is measurable with respect to 
V (g) B(R) <8> B{R d ) and the terminal condition £ is simply a i? m -valued .^-measurable 
random variable, k is defined on A x A and scalar valued. 

By a solution to RBSDE (1.1) we mean a triple (Y,Z,K) = {Y(t), Z(t), K(t)} te[0:T] 
of predictable processes with values in R m x R mxd x R m such that P-a.s., t -» Y(t) 
and i — > if(i) are continuous, t — ► Z(t) belongs to L 2 (0, T), t — > ^(i, Y^t), Zi(t), i) 
belongs to L l {d,T) and P-a.s., RBSDE (1.1) holds. 

S 2 (R m ) or simply S 12 denotes the set of P m -valued, adapted and cadlag processes 
{ Y (t)}te{o,T\ such that 

||Y|| 52 :=E[ sup |y(t)| 2 ] 1/2 <+«>. 
te[o,T] 

(S 2 , || • H52) is a Banach space. 

M 2 (R mxd ) or simply M 2 denotes the set of (equivalent classes of) predictable 
processes {^(t)}te[o,T] with values in R mxd such that 

\\Z\\ M 2 := E[ [ T {Zis^ds] 1 / 2 < +00. 
J 

M 2 is a Banach space endowed with this norm. 

N 2 (R m ) : = {K = (Ki, • ■ ■ , if m ) T € S 2 : for any i € A, *Q(0) = 0, 
and t — > -ftTj(t) is increasing }, 

where T means transpose. (iV 2 , || ■ H52) is a Banach space. 

2.2 Existence 

In this subsection, we prove the existence result for RBSDE (1.1). We assume the 
following Lipschiz condition on the generator. 

Hypothesis 2.1. (i) ^(-,0,0) := (^(-,0,0, 1), ■ ■ ■ , tp(-, 0, 0, m)) T belongs to M 2 . 

(ii) There exists a constant C > 0, such that, P-a.s. for each (t,y,y',z,z',i) € 
[0, T) x R x R x R d x R d x A, 

IV-fty,*,*) <C(|y-y'l + k-^l)- 
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We make the following assumption on k which is standard in the literature of 
optimal switching. 

Hypothesis 2.2. (i) For any G A x A, k(i,j) > 0. 
(ii) For any G A x A x A, 

k(i,j)+k(j,l) > k(i,l). 

We are now in position to state the existence result. 

Theorem 2.1. Let the Hypotheses 2.1 and 2.2 hold. Assume that £ G L 2 (£l, T T , P; R m ) 
takes values in Q. Then RBSDE (1.1) has a solution (Y, Z,K) in S 2 x M 2 x N 2 . 

We first sketch our proof. 

Sketch of the Proof: The proof is divided to four steps. In Step 1, we introduce 
the penalized BSDEs whose existence and uniqueness follows from the classical result. 
In Step 2, we state some (uniform) a priori estimates for the solutions of penalized 
BSDEs, whose proof will be given in the next subsection. In Step 3, we prove the 
(monotone) convergence of these solutions. Finally, in Step 4, we check out the minimal 
boundary condition. 

Proof of Theorem 2.1 : Step 1. The penalized BSDEs. 

For any nonnegative integer n, let us introduce the following penalized BSDE: 

Y?(t) = & + j\{s,Y l n {s),Z^{s) 1 t)ds 

-nJ2 / (Y^s)-ms)-k(i,l)) + ds (2.1) 
l=i Jt 

- Z?{s) dW(s), t G [0, T], i G A. 

Note that when I = i, we have, in view of Hypothesis 2.2 (i), 

(Y?(s)-Y?(s)-k(i,l)) + = 0. (2.2) 

From the classical result of Pardoux and Peng [14], for any n, BSDE (2.1) has a 
unique solution (Y n , Z n ) in the space S 2 x M 2 . 

Step 2. A priori estimates. 

The following lemma will play a crucial rule in the proof of Theorem 2.1. 

Lemma 2.1. Let the Hypotheses 2.1 and 2.2 hold. Let us also assume that £ G 
L 2 (Q, Tt-i P-, R m ) takes values in Q. Then there exists a constant C > (indepen- 
dent of n ), such that 

\\Y n \\ S 2 + \\Z n \\ M 2<C (2.3) 

and 

n 2 E f ((Y?(s) - Y?( 8 ) - k(i,j)) + f ds < C. (2.4) 
J o 
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However, the proof of this lemma is quite lengthy and delicate. We relegate it to 
the next subsection. 

Step 3. Convergence of solutions {Y n , Z n } of penalized BSDEs. 
First, for each n, we introduce a function ip n as follows: 

m 

^ n (t,y,z,i) -WMZirf-n^iyi-yt-kfrl))*, (t,y,z,i) G [0, T] x R m x R mxd x A. 

1=1 

Since ip n (-, z,i) depends only on Z{ and 

4> n {t, y + y', Zi,i) > ip n (t, y', Zi,i) 
for any y G R m such that y > and the i — th component j/j = 0, it is easy to check 

m 

-4(y- , V n (t, y + + z) - r + \t, y', z')) < 2 £ Xj/i <o|^ - ^'| 2 + ^V| 2 , P - a.*., 

i=l 

for a constant C > 0. We note that this is the inequality (5) in [12]. Applying the 
comparison theorem for multi-dimensional BSDEs (see Hu and Peng [12, Theorem 
2.1]), we deduce that for any nonnegative integer n, 

Y?(t) > y™ +1 (t),Vi <E A,t £ [0,T]. (2.5) 

For a.e. t and P-a.s. u, {Y n (t,uj)} n admits a limit, denoted by Y(t,u>). Moreover 
from the a priori estimate (2.3) and Fatou's lemma, we have 

sup E\Y(t)\ 2 < C. (2.6) 

te[o,T] 

In view of the fact that Y^t) < Y t n (t) < Y?(t) with i € A, applying Lebesgue's 
dominated convergence theorem, we have 



lim E f \Y n {s)-Y{s)\ 2 ds = 0. 

n->oo J Q 



(2.7) 



Now we prove that {(Y n , Z n )} n is a Cauchy sequence in the space S 2 x M 2 . For 
this purpose, we apply Ito's formula to \Y™(t) — Yf(t)\ 2 to obtain 

\ Y -(t) - Y[{t)\ 2 + £ \Z?(s) - Zf(s)\ 2 ds 

= 2 [(Yr{s) - Y?(s)ms, Y?{8), Z?(s),i) - Y,(s, Yf(s), Zf(s),i)) ds 

-2 J\y?(s) - Y*{8))n{Y?{8) - Y?{s) - k(i,j)) + ds 

+2^V«( S ) " Yf(s))p(Yf(s) - Yf(s) - k(i,j)) + ds 

-2^V«( S ) " Y?{s)){Z?{s) - Zf(s))dW(s),i G A. (2.8) 
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Putting t = and taking expectation in the last equality, we get for i € A, 

E\Y?{0) - Yf (0)| 2 + E f \Z?{a) - Zf(s)\ 2 ds 

Jo 

= 2E l T {Y l n {s)-Yl{s))^{ S ,Y l n { S ),Z^s)^)-^{s,Yl{s),Zf{ S ),t))ds 
Jo 

-2E f(Yt(s) - Y?{8))n<Y?{ 8 ) - Y?{s) - k^j)) + ds 
Jo 

+2Ej\Y l n ( S ) - Yf(s))p(Yf(s) - Yf(s) - k(i,j)) + ds 

< CEj*(Y?{8) - Y[(s)fds + ±e£ \Z?{s) - Zf(s)\ 2 ds 

T x I / rT 2 x , 



+ \Ej \Y™(s)-Yf(s)\ 2 ds) \Ej n 2 {{Y l n {s)-Yp{s)-k^j))+Yds 

+ (e£ \Y?{s) - Yf(s)\ 2 ds^j 2 [eJ\ 2 (p?{8) - Yj{s) - k(z,j)) + Yds 
From (2.4) and (2.7), we have 

lim E [ T \Zl l (s)-Zf(s)\ 2 ds = 0,i£ A. (2.9) 

n,p->oo J 



Now, we define the increasing process Kf as follows: 

r t m 



K{t):=n y2(Y l n (s)-Y l n (s)-k(i,l)) + ds, te[0,T],ieA. (2.10) 
Jo i=i 

From the penalized BSDE (2.1), we have 

K?(t) = Y?(t)-Y?(0)+ fl,(s,Y?{8),Z?(a),i)d8- f Z? (s)dW (s) , i e A. (2.11) 

Jo Jo 

We denote by Z the limit of Z n in M 2 . Set 

Ki(t):=Y(t)-Y t (0)+ [ i/>(s,Yi(s),Zi(s),i)ds- I Z t (s)dW{s),i e A. (2.12) 
We have 

lim ||Z n -Z|| M2 =0. 

n^oo 

Going back again to (2.8), we deduce that for i E A, 
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E sup \Y?{t)-Y?{t)f 

0<t<T 



< 2E f \Y^s)-Yl( S )m S ,Yr( S ),Z^sU)-^{s,Y[{s),Zf{s)Md S 
Jo 

+2E [ T \Y?(8) - Y?(s)\n(Y?( 8 ) - Y?( 8 ) - k(z,j)) + ds 
Jo 

+2E jf \Y;\s) - Yf(s)\p(Yf(s) - Yf{s) - k(i,j)) + ds 

+2E sup | [ T (Y?(s)-Yf(s))(Z?(s)-Zf(s))dW(s)\. (2.13) 

0<t<T Jt 

The last term of the right hand side of the last inequality is less than or equal to 
the following quantity: 



CE QT \(Y?{s) - Y?{s)){Z?{s) - Zf(s))\ 2 ds^j 2 

< CE[ sup \Y-(t)-Yf(t)\ [ T \Z?(s)-Z?(s)\ 2 ds) 

\0<t<T Jo J 

< \E sup \Y?(t)-Y?(t)\ 2 + CE [ T \Z?( S )-Z!(s)\ 2 ds. (2.14) 

^ 0<t<T Jo 

Combining (2.13) and (2.14) and taking into consideration (2.7) and (2.9), we de- 
duce that {Y n } n is a Cauchy sequence in S 2 , which means in particular that 

lim \\Y n - Y\\ S 2 = 0. 

n— >oo 

Consequently, Y is a continuous process. 

From (2.11), (2.12), and the following fact that 

\\Y n -Y\\s» + \\Z n - Z\\ M 2 -> 0, asra^oo, 

we deduce immediately that 

lim \\K n - K\\ S 2 = 0. 

Hence, K £ N 2 , and (Y, Z,K) satisfies the first relation in RBSDE (1.1). 
Finally, from the a priori estimate (2.4), we have 

E £ {(Y?(s) - Y?(s) - k(i,j)) + ) 2 ds < ^, i,j € A. 
Jo n 

Letting n — > oo, we deduce 

E / ((Yi(s) - Yj(s) - k(i,j)) + f ds = 0, i,j e A. 
Jo 



Hence, 

YitfKYjW + kiiJ), s G [0, T],i,j G A, (2.15) 
which are equivalent to the following: 

P-a.s. Y(s)€Q, VsG[0,T]. 

Step 4. The minimal boundary condition. 

Let us first state the following lemma whose proof is at the end of this subsection. 

Lemma 2.2. Let the Hypotheses 2.1 and 2.2 hold. Let us also assume that £ G 
L 2 (S7, Tt-, P; R m ) takes values in Q. We have, for any integer n, 

J o (w(s)-mm\Y?{s) + k{i,j)]) dK?(s) = 0, i G A. (2.16) 

Now, we can take the limit in (2.16) by letting n tend to +oo and applying Lemma 
5.8 in [10] to get the following 



a 



Yi(s) - min [Yj(s) + k(i,j)U dK^s) =0, % G A, (2.17) 



which, together with (2.15), yields the minimal boundary conditions. 

The proof of Theorem 2.1 is now complete. □ 
Proof of Lemma 2.2 For ieA, the left hand side of the last equality is equal to 

the following sum 



m „t 

/ mm{(Yr(s)-Y J n (s)-k(z,j)y(Yr(s)-Yr(s)-k( l ,l)) + }ds. (2.18) 



We claim that the integrand of the l-ih integral is equal to zero for / G A. 
In fact, it is immediate for the case of / = i. For the case of Z ^ i, the integrand is 
the minimum of the following m — 1 nonnegative quantities: 

" Y?(s) - k(i,j))- (Y?(s) - Y?(s) - k(i, l)) + , j G A and j + i, (2.19) 

whose l-th term is zero due to the fact that 

(Y?(s) - Y?(s) - fc(i,Z))- (Y^s) - Y?{a) - k(i,l)) + = 0, (2.20) 

and therefore, it is zero. The desired result then follows. □ 

2.3 Proof of Lemma 2.1 

In this subsection, we prove Lemma 2.1. 
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For i,j G A, applying Tanaka's formula (see, e.g. [17]) to (Y^ l (t) — YJ l (t) — k(i,j)) + , 
we have 

m „r 

(Y?(t)-Y?(t)-k(i,j)) + +nJ2 / X £+ (*)(!?(*) -^-MM))^ 

~ n J2 Xc+ n (s)(Y?(s)-Yr(s)-k(j,l))+ds + - dL^(s) 
1=1 Jt l3 ' n Jt 

= I X C+ {s)^{s,Y?{s),Z?{s),i)-^{s,Y?{s),Z^s),j))ds 
- f Xc+ (s)(Z?(s)-Z?(s))dW(s), 

(2.21) 

where for i,j £ A, 

4,„ == iM : Y?(s) > Y?(s) + k(i,j)}, (2.22) 

and L™- is the local time of the process Y t n — YJ 1 — k(i,j) at 0. 

Applying Ito's formula to ^(Y^ n (i) — Yp(t) — k(i,j)) + ^j and taking into consider- 



ation 



j\Y l n {s)-Y^{s) - k{t,j)) + dL^{s) = 0,V/ G [0,71, (2.23) 



we have 

{(Y:?(t)-Y?{t)-k{i,j))+) 2 

r-T 



+2nJ t ((Y;\ S )-Y?{s)-k(t,j)) + f ds 

+ f Xc+ (s)\Z?{s)-Z](s)\ 2 ds 

Jt l] ' n 

= 2j\Y t n (s)-Yj;(s)-k( i ,3)) + ty(s,Y?( 8 ),Z?(s),i)-4(s,Y?(s),Z?(s),j)] ds 
-2 f(Yr(s) -Y?(s) -k(i,j))+(Z?(s) -Z?(s))dW(s) 

Jt rp 

+2nj t (Yl>(s) - Y?{s) - k{i,j)) + {Y?{s) - Y?(s) - k(j,i)) + ds 
+2n £ [ T (Y;\s)-Yj\s)-k( i ,j)) + 

x [<y?{8) - ma) - hj, /))+ - (xns) - m S ) - Hi, /))+] ds. 

(2.24) 

We claim that the integrands of the integrals in the last two terms of (2.24) are all 
less than or equal to zero. In fact, since 

{(yi, • • • , y m ) T €R m :yi- yj - k{i,j) > 0, Vj - Vi - k{j, i) > 0} = 

due to the fact that 

k(i,j) + k(j,i) > 0, 
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we have 



(Y?( 8 ) - Y?( 8 ) - k{i,j)) + {Y?{s) - Y?{s) - k{j,i))+ = 0, i,j e A. (2.25) 

On the other hand, for l,i,j G A, taking into consideration both Hypothesis 2.2 
(ii), i.e., 

k(i,j)+k(j,l) > k(i,l), 

and the elementary inequality that — x^ < [x\ — X2) + for any two real numbers x\ 
and X2, we have 

(17(a) - Y?(s) - k(i,j))+ [(Y?(s) - Yp(s) - k(j, l))+ - - Y?(s) - k(i, 1))+} 

< (Y?(s) - Y?(s) - k(i,j)) + {Y?{s) - Y t n (s) - k(j, I) + k(i, l))+ 

< (Y?(s) - Yf(s) - k(i, I) + k(j, l))+(Yp(s) - *700 - Hj, I) + k(i, /))+. 

(2.26) 

The last term of the last inequality is zero, since 

{(yi, ■ • • , Vmf eR m : yi - yj - k(i, I) + k(j, I) > 0, Vj - Vi - k(j, I) + k(i, I) > 0} = 0. 
Concluding the above, we have 
E((Y?(t)-Y?(t)-k(i,j))+) 2 

+2nE£ ((Y t n (s) - Y?{s) - k{t,j)) + f ds 

+E [ T Xc+ (s)\Z?(s)-Z?(s)\ 2 ds 

Jt lJ ' n 

< 2E j 1 \y?{ 8 ) -Y?{s) - k{i,j)) + W(s,Y?{s),Z?(s),i) - ^s,Y J n (s),Z^s),j)\ds. 

(2.27) 

In view of Hypothesis 2.1 on the function tp, we have 

\^(s, Y?(s), Z?(s), i) - ^(s, Y?{s), Z?(8),j)\ 

< \^(s, 57(a), Z?(s), i) - </>(a, 17(a), Z?(s),j)\ 
+ \iP(s,Yi>(s),Z?(s),j)-^(s,Y?(s),Z?(s),j)\ 

< 0(^(3, 0,0)| + |Y/'(s)| + \Z?{s)\ + \Y?(s) - YJ\s)\ + \Z?(s) - Z?(s)\) 

< C(l + ms,0,0)\ + \Yi>(s)\ + \Z?(s)\ 

+ \Y?(s) - Y?( 9 ) - k(i,j)\ + \Z?(s) - Zf(s)\). 

(2.28) 
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Consequently, we have 

E{(Y?{t)-Y?{t)-k(i,j)) + ) 2 

cT 



+2nEj t ((Y?(s) - Y?(s) - k{i,j)) + ) 2 ds 

+E f Xc+ (s)\Z?(s)-Z^(s)\ 2 ds 

Jt %3 ' n 

< CE £ \(Y?(8) - Y?{s) - k(z,j)) + \ 2 ds 



t ^ n J (2.29) 

T 



it 



+\ E I x c ± ( s )(i + IV( S ,o,o)| 2 + |y™( s )| 2 + |zr(,)| 2 

- Y?(s) - k(i,j))+\ 2 + \Z?(a) - Z?(s)\ 2 ) ds. 
Applying Gronwall's inequality we deduce easily that 

E ((Y{\t) - Y?{t) - k{i,j))+) 2 < C (l + E f\ ct Js) (|Y«(s)| 2 + \Z?(s)\ 2 ) ds^j , 

nE f ((Y?(s)-Y?(s)-k(i,j)) + ) 2 ds + E [ T Xc+ {s)\Z?{s)-Z?{s)\ 2 ds 
Jo Jo %3 ' n 

<c(\ + E £ Xc+Js)[\Y?{8)\ 2 + \Zf{s)\ 2 ] ds^j . 

(2.30) 

Going back to (2.24) and applying Burkholder-Davis-Gundy's inequality we obtain 
E[ sup ((Y-(t) - Y?{t) - fc(i,i))+) 2 ] <c(l + E f Xc± (s)[\Y^(s)\ 2 + \Z?(s)\ 2 ] ds 

0<t<T V Jo t3 ' n 

(2.31) 

On the other hand, from (2.27), we deduce that, 

E ((YT(t) - Y?(t) - k(i,j))+) 2 + 2nE £ ((Y^s) - Y?{s) - k(i,j)) + ) 2 ds 
+E [ Xc+ (s)\Z?(s)-Z?(s)\ 2 d S 
< (n + C)E £ ((Y?(a) - Y?{s) - k(i,j)) + ) 2 ds 

+^-E ! T Xc ± W(l + lV>M,0)j 2 + \Y?(s)\ 2 + \Zf{s))\ 2 + \Z?(s) - Z-{s)\ 2 )ds. 



It 

This shows that, for sufficiently large n 

cT 



n 2 E jT ( ( y« (s ) _ Y J\s) - k(i,j)) + ) 2 ds<c(l + E jf [|y«(,)| 2 + \Zf(s)\ 2 ] ds^j . 



(2.32) 
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Finally, applying Ito's formula to \Y™{t)\ , we obtain: 

pT 



\Y{\t)f + J t \Z?{s)\ 2 ds 
-2 J Zi(s)dW(s). 



1>(8,Y?(8),Z?(8),i) -J2n(Y?(s) - Y t n (8) - k(i,l)Y 



By taking expectation and using the elementary inequality: 

2ab < -a 2 + eb 2 , Ve > 0, 
e 

we deduce that, for any e > 0, 

E^^ + E £ \Z?{s)\ 2 ds 

< EM 2 + 2E £ \Y?( S )\- [Ms, Y?{s), Z?{s),i)\ 

m 

+ Y^n(Y^(s)-Y l n ( S )-k( i ,l)) + ] ds 
l=i 

< c + 2E / T |y«( s )|.[c(|^,o,o)| + |y«( s )| + |zr( s )|) 



+ J2n(Y?(s)-Y l n (8)-k(i,l)) 
l=i 

C + C e E [ \Y?(s)\ 2 ds + eE f \Z?(s)\ 2 ds 
Jt Jt 

rT m 

+eE / n 2 J2{(Y^)-yi n (s)-Ki,l)) + ) 2 ds 
Jt i=i 

f \Yr{s)\ 2 ds + CeE [ T \Z?(s)\ 2 ds. 
Jt Jt 



< C e + C £ E 



ds 



(2.33) 



Here, C e > denotes a constant which depends on e and may vary from line to line. 
Therefore, 



E\Y?(t)\ 2 + E^ \Z?(s)\ 2 ds<C. 
From (2.32), we obtain (2.4), and from (2.33), we deduce 

\\Y n \\ s i<C. 
The proof of Lemma 2.1 is now complete. 



(2.34) 

(2.35) 
□ 
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3 Uniqueness 



In this section, we prove the uniqueness by a verification method. Let (Y,Z,K) be a 
solution in the space (S 2 ,M 2 ,N 2 ) to RBSDE (1.1). We will prove that Y is in fact 
the (vector) value for an optimal switching problem of BSDEs. For this purpose, we 
introduce the following optimal switching problem. 

Let {^jj-j^o ^ e an increasing sequence of stopping times with values in [0, T] and 
Vj, aj is an J-g. -measurable random variable with values in A, and x is the indicator 
function. We assume moreover that for P-a.s. uj, there exists an integer N(lu) such 
that 6 N = T. 

Then we define the admissible switching strategy as follows: 

N 

a(s) = a oX {e }(s) +^ j a j - lX (e j _ 1 fi j ]{s)- (3.1) 

3=1 

We denote by A the set of all these admissible switching strategies and by A 1 the 
subset of A consisting of admissible switching strategies starting from the mode i. In 
the same way, we denote by At the set of all the admissible strategies starting at the 
time t (or equivalently 6q = t ) and by A\ the subset of At consisting of admissible 
switching strategies starting at time t from the mode i. 

For any a(-) G At, we define the associated (cost) process A a ^ as follows: 

N-l 

A<\s) = fc(ai-i,a;)Xfo,n(s), s G M- ( 3 - 2 ) 
j'=i 

Obviously, A a( --\-) is a cadlag process. 

Now we are in position to introduce the switched BSDE: 

U(s) =Z a{T) +A a ^(T)-A<\s)+ f \fj(r,U(r),V(r),a(r))dr- f V(r)dW(r),s G [t,T\. 

J S J S 

(3.3) 

This is a (slightly) generalized BSDE: it is equivalent to the following standard BSDE: 
17(a) = £ a(T) +A a( --\T)+ I i>(r,U{r)-A a( --\r),V(r),a(r))dr- ! V(r)dW(r),s G [t, T] 

J S J S 

(3.4) 

via the simple change of variable: 

U(s) = U(s) + A<\s), V(s) = V(s), s£[t,T\. 

Hence, BSDE (3.3) has a solution in S 2 xM 2 . We denote this solution by (U a ^,V a ^). 
Note that U is only a cadlag process. 

The optimal switching problem with the initial mode i G A is to minimize U a ^ (t) 
subject to a(-) G A\. 

The assumptions required for the uniqueness will be slightly stronger than those 
needed for existence. We keep the same assumption on ip and we assume the following 
for k. 
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Hypothesis 3.1. (i) For any £Ax A, k(i,j) > 0. 

(ii) For any I) £ A x A x A such that i ^ j and j ^ I, 

k(i,j) + k(j,l) > k(i,l). 

We have the following representation for the first component of the adapted solution 
to RBSDE (1.1), which immediately implies the uniqueness of the adapted solution to 
RBSDE (1.1). 

Theorem 3.1. Let us suppose that the Hypotheses 2.1 and 3.1 hold. Let us also 
assume that £ € L 2 (Q,J 7 T,P',R m ) takes values in Q. Let (Y,Z,K) be a solution in 
{S 2 ,M 2 ,K 2 ) to RBSDE (1.1). Then 

(i) For any o(-) G A\, we have: 

Yi{t)<U<\t), P-a.s. (3.5) 

(ii) Set 0q = t, ctQ = i. We define the sequence {6*,a*}°? =1 in an inductive way as 
follows: 

0* := M{s > 9*_ x : Y aU (s) = min {Y t (s) + I)} A T, (3.6) 

and a* is the To* -measurable random variable such that 

Y aU (9*) = Y a *{6*) + fc(«*-i,«*), 

with j = 1, 2, • • ■ . 

Then, P-a.s. uj, there exists an integer N(lu) such that 6* N = T . And the following 
switching strategy: 

N 

a*{s) = iX{t}{s) +^2 a j-iX(9*_ 1: e*](s), (3.7) 
j'=i 

is admissible, i.e., a*(-) € A\. Moreover, 

Y(t) = u°'U(t). 

(Hi) We have the following representation for Y{t): 

Yi(t) = essinf U a ^{t), i e A,t € [0,T]. 
a(-)eA\ 

RBSDE (1.1) has a unique solution. 

Proof. Without loss of generality, we will prove (i) and (ii) for the case of t = 0. 
Otherwise, it suffices to consider the admissible switching strategies starting at time t. 
(i) We define 

N 

Y<\s) = Y,^-^m-iA)^) + L(T)X { T}(s), (3.8) 
i=i 

N 

Z<\s) = YsK-Mx^fi^s), (3.9) 

N i-B-As 

K<)(s) = £ / dK ai _,{r). (3.10) 

i=l JSi-ihs 
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Noting that Y a ^(-) is a cadlag process with jump Y ai (9i) — Y ail (6i) at 0j, i = 
1, • • ■ , AT — 1, we deduce that 

y«(-)( a )_ y«(0(o) 

= E / hV'(r,^ i _ 1 (r),^ i _ 1 (r),a i _ 1 )dr + ^ i _ 1 (r)dW(r) + ^ ai _ 1 (r)] 

i=l •M-lAs 
N-l 

i=l 

= [\-ip{r, Y a ^{r), Z<\r), a{r))dr + Z^rJdW^r) + dK a «(r)] 
JO 

where 

a<\ s ) = + H^-i,ai) - Y^e^x^nis), (3.11) 

1=1 

and it is an increasing process due to the fact that 

Y(t)eQ, VtG[0,T]. 
Consequently, we conclude that (Y a ('\ Z a ^) is a solution of the following BSDE: 

= L(T) + A a{ - ] (T) - A a ^ (s) - [(K a ^ (T) + A a ^ (T)) - (K a ^ (s) + A a U ( s ))] 
+ f i){r,Y<\r),Z a{ --\r),a{r))dr - f Z<\r)dW{r), s G [0, T}. (3.12) 

J s J s 

Since both K a ^ and A a ^ are increasing cadlag processes, from the comparison 
theorem, we conclude that 

y-«(0(o) < ?7 a (-)(o), 

which implies that 

Yi(0) < U a ^{0). 

(ii) Let us first claim that if < 9\ < Q\ < T, then there exists a constant c > 
such that 

\Y(9* 2 )-Y(ei)\>c. 
To prove this claim, we introduce the following subsets of Q: for i / j, 

Bij := {(yi, • • • ,y m ) T eR m : yi = Vj + k(i,j)} n Q. 

We assert that for i ^ j and j / Z, n Bj i = 0. 

In fact, if there exists an element (yi, • • ■ , y m ) T G -Bjj n S^;, then 

y; = yj + fe(i, j) and yj = y x + /). 
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We deduce then 

Vi = yi + k(i,j) + k(j,l) > yi + k(i,l), 

which contradicts the fact that (yi, • ■ ■ , y m ) T G Q. 

Hence, the distance between Bij and Bji is strictly positive, 

dist(Bij,B jt i) > 0. 

We set 

c := min dist(Bij, Bj i) > 0. 
We return to the proof of the claim. From the definition of (0*, a*) and (6*2, a^), 

€B ita . fmdY(eZ)€B altt4 , 

which implies that 

\Y(6* 2 )-Y(6l)\>c, 

and the proof of the claim is finished. 

In the same way, if 0* < 0% < ■ ■ ■ < 6»*_ x < 9* < T, then 

|y(e;)-n<9;_i)l>c. 

On the other hand, as Y satisfies (1.1), it is easy to check that 



E 



£ \Y{6*) - Y{6U) 

3=1 



< oo. 



As a consequence, there exists N(u) such that 0^ = T. 
Finally, from the choice of a*(-), 

K a *V+A a *U = 0. 

We conclude from (3.12) that 

Y a *(-)( ) = [T*(-)(0), 

which implies that 

y.(0) = !/«*(■) (0). 

(iii) The representation for Y is a combination of both assertions (i) and (ii) . This 
gives the uniqueness of the first component of the adapted solution, and the uniqueness 
of the other two components of the adapted solution follows then. □ 
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4 Optimal switching of functional SDEs 



In this section, we study the optimal switching problem. In order to ensure the exis- 
tence of optimal switching strategy, we use the weak formulation of the problem. Let 
(Q, T, P) be a complete probability space and let {Tt,t > 0} be a filtration satisfying 
the usual conditions. The process W is an {J^,t > 0}-Brownian motion on R d defined 
on (Q,F,P). 

Consider the switched equation 

X<\t) = x + f a(s,X a ^)[dW{s) + b(s,X a( -\a{s))ds], t e [0,T] (4.1) 
Jo 

and the cost functional 



J(a(-)) = E 



T 



l{s,X a( -\a{s))ds 
o 



+ E 



7V-1 



^ k(oti-i,ai) 



i=l 



(4.2) 



The switching problem is to minimize the cost J(a(-)) with respect to a(-), subject to 
the state equation (4.1). 

In the above, xq is a fixed point in R d . a, b and I are defined on [0, T] x C([0, T];R d ), 
[0,T] x C([0,T];R d ) x A and [0,T] x C([0,T];R d ) x A, respectively, with values in 
R dxd , R d and R, respectively. As in Section 3, a(-) is an admissible (J-i)t>o-adapted 
switching strategy, A is the set of all the admissible {J-~t,t > 0}-adapted switching 
strategies and A 1 is the subset of A consisting of the admissible {Tt,t > 0}-adapted 
switching strategies starting from the mode i. We assume that a, b(-,i) and are 
progressively measurable functionals on C([0, T];R d ) in the following sense: 

Definition 4.1. Let C([0, T]; R d ) be the space of continuous functions x : [0, T] — > 
For < t < T, define Qt ■= cr(x(s) : < s < t), and set Q := Qt- A progressively 
measurable functional on C([0,T]; R d ) is a mapping jjl : [0,T] x C([0,T]; R d ) —> H ( 
H is some Euclidean space) such that for each fixed t G [0, T], fi restricted to [0,t] x 
C([0,T];R d ) is B([0,t]) ® Q t /B{H)-measurable. 

We assume that k satisfies Hypothesis 3.1. And we assume also that a, b and I 
satisfy the following hypothesis. 

Hypothesis 4.1. (i) a, b(-,-,i) and l(-,-,i),i € A, are progressively measurable func- 
tionals on C([0,T];R d ). 

(ii) There exists a constant (3 > such that V(t,x,x',i) £ [0, T] x R d x R d x A, 

\b(t,x,i) - b(t,x',i)\ + \a(t,x) - a{t,x')\ + \l(t,x,i) - l(t,x',i)\ < 0\\x - x'\\ c{[ ^ t] . Rd y 

(Hi) a has a bounded inverse, 
(iv) b is bounded. 

Let (Y, Z, K) be the unique solution in (S 2 , M 2 ,N 2 ) of the following RBSDE: 

f Yi(t) = f 1>(s,X,Zi(s),i)ds- [ dKi(s)- f Zi{s)dW{s), 
Jt Jt Jt 

Yi(s) < mm{Y j (s) + k(i,j)}, (4.3) 

J (Y l (s)-mm{Y J (s)+k(i,j)}^dK l (s)=0, i G A, 
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where V is defined as follows: V(t,x,z,i) G [0, T] x C{[0,T];R d ) x R d x A, 

^(t, x, z, z) := l(t, x, i) + (z, b(t, x, i)), 
and X is the solution to the following functional SDE: 

X(t) = x + f a(s,X)dW{s), t€[0,T\. (4.4) 
Jo 

Theorem 4.1. Let the Hypotheses 3.1 and 4-1 hold. Then 

(i) For any a(-) G .A*, we have: 

J(a(-)) > Yi(0). (4.5) 

(ii) There exists an optimal switching strategy a* , and a weak solution (P*, W*,X*), 
such that 

X*(t) = x + [* a(s,X*)[dW*(s)+b(s,X*,a*(s,X*))], t e [0,T], (4.6) 
J o 

and 

J(a*(-))=Y t (0). 
Proof, (i) For any a(-) G A 1 , we set 

dP:=expt- J b(s,X a (-\a(s))dW(s) \b(s, X<\ a(s))\ 2 ds^ dP. 

Then W>(t) = W(t) + J *6(s,X a (-),a(s))ds is a Brownian motion under the new prob- 
ability measure P. Let (Y, Z, K) be the solution of the following RBSDE: 

Yi(t) = f 4>(s,X<\z l (s),i)ds- [ dKi(s)- f Zi( 8 )dW(s), 

Jt Jt Jt 

Yi{s) < mm{Y j (s) + k(i,j)}, (4.7) 

J (Y i (s)-mm{Y j (s) + k(iJ)}^dK i (s) = 0, i G A. 

Note that since X a ^ solves (4.1), we have 

X<\t)=x + f a{s,X a ^)dW(s), t€[0,T\. 
Jo 

By a classical argument of Yamada-Watanabe, for RBSDE (4.3), the pathwise unique- 
ness implies the uniqueness in the sense of probability law (see, e.g. [5], for a proof in 
the framework of BSDE). Hence, we have 

Yi(0) = Yi(0), ieA. (4.8) 
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Recalling the cost process A a ^ denned by (3.2), and denning Y<\ Z a ^ , K a ^ 
and A a ^ in the same way as (3.8), (3.9), (3.10) and (3.11), we deduce in the same 
manner as in Section 3 that (Y a U , Z a ^) is the unique solution of the following BSDE: 

= A<\T)-A<\t) - [{K<\T) + A<\T)) - {K<\t) + A<\t))] 

+ / ij(s,X a (-\z a( --\s),a(s))ds- [ Z a ^(s)dW(s), t € [0,T\. (4.9) 
Jt Jt 

Since, both K and A are increasing, we deduce that 
y«(0(o) 

rT _ r T _ 

< A<\T)+ iP{s,X<\z<\s),a{s))ds- Z a ^(s)[dW(s) + b(s,X<\a(s))} 
Jo Jo 

= A<\T)+ [ l{s,X<\a{s))ds- f T Z<\s)dW{s). (4.10) 
Jo Jo 

From the definition of Y a ('\ 

Yi(0) = y a ^(0). (4.11) 

From (4.8), (4.11) and (4.10), 

y(0) < A<-\T)+ [ T l(s,X a( --\a(s))ds- f T Z<\s)dW{s). 
Jo Jo 

Taking expectation with respect to P, we have 

y(0) < J(a(-)). 

(ii) Let X be the solution of SDE (4.4) and (Y, Z, K) be the solution of RBSDE 
(4.3). Then, Y is adapted to the filtration = , due to (4.4) and Hypothesis 4.1 
(iii). 

Set 0q = 0, ctg = i- We define the sequence {6*,a*}'jL 1 in an inductive way as 
follows: 

6* := inf {s > 9*^ : Y a * (a) = min {Y,(s) + fe(aj_ 1 , 0} A T, (4.12) 

and a* is the J^* -measurable random variable such that 

Y aU (0*)=Y a *(9*) + k(a*_ 1 ,a* j ), 

with j = 1,2, ••■ . Then, P-a.s. u>, there exists an integer N(u) such that 6* N = T. 
And we define the switching strategy a* as follows: 

N 

a*(s) := iX{o}(s) + J2uj-iX(e*_ 1 ,e*](s), (4.13) 
i=i 

is admissible, i.e., a*(-) £ A 1 , a* is adapted to the filtration = , i.e., a*(t) = 
a*(i,X). Setting 

(LP* :=expQf b(8,X,a*(8,X))dW(s) \b(s,X,a*(s,X))\ 2 ds^ dP, 
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then 



W*{t) := W(t) - [ b{s,X,a*(s,X))dW{s), t€[0,T\, 
Jo 



is a Brownian motion under the probability measure P*, and (P* ,W* , X) is a weak 
solution of (4.1). 

Computing Y a *(') as in Section 3, we deduce that 

y«*(-)( ) 

= A a * ( -\T) - {K a '^{T) + 2 a *(-)(T)) 

+ f if>(s,X,Z a *(-Xs),a*(s))ds - [ Z a *^(s)dW(s), (4.14) 
Jo Jo 



where 



JV-1 



A a*(-) 
Z a*(.) 



3=1 
N 

zZ Ya *-M X i e t-i^t)^ + Y a*(T)(s)X{T}(s), 
i=l 
N 



n r e*As 



N-l 



i=i 



From the definition of a*(-), K a *(') = and ^4 a *(') = 0. Hence, from (4.14), it follows 
that 

y«'(0(o) 

= ^ a *»(T)+ / Z(s,X,a*(s))ds- / Z a *»(s)(W*(s). (4.15) 

./0 JO 

Taking the expectation with respect to P*, we conclude the proof. □ 

5 System of variational inequalities 

In this section, we will show that the RBSDE studied in Sections 2 and 3 allows us to 
give a probabilistic representation of the solution to a system of variational inequalities. 
For this purpose, we will put RBSDE (1.1) in a Markovian framework. 

Let b : [0,T] x R d -> R d and a : [0,T] x R d -> R dxd be continuous mappings. We 
assume: 
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Hypothesis 5.1. There exists a constant C > 0, such that for all t G [0, T], and 

(x,x') G R d x R d , 

\b(t,0)\ + \a(t,0)\<C, 



and 



\b(t,x) - b(t,x')\ + \a(t,x) - a(t,x')\ < C\x - x'\. 



For each (t,x) G [0, T] x R d , let {X^ x {s);t < s < T} be the unique .Revalued 
solution of the SDE: 

X t ' x (s)=x + f S b(r,X t ' x (r))dr+ f a{r, X l ' x {r))dW{r), s G [t,T]. 
Jt Jt 

We suppose now that the data £ and tp of RBSDE (1.1) take the following form: 

& = g(X^(T),i), 
i>{s,V,z,i) = ^{s,X^ x (s),y,z,i), 

where g and ■0 are given as follows. 

Hypothesis 5.2. (i) For each i G A, the function g(-,i) £ C(R d ) and has at most 
polynomial growth at infinity, 
(ii) For each i £ A, 

ip{-,i) : [0,T] x R d x R x R d ^ R 

is jointly continuous and there exist two constants C > and p > such that 

|V»(t,x,0,0,i)| < C(l + 

|V>(t,x,y,z,i) - ip(t,x,y',z',i)\ < C{\y - y'\ + \z - z'\), 
for t G [0,T],x,z,z' G R d ,y,y' G i? ; and i G A. 



(Hi) Vx G 



g( x ) : = (a(x, !),-•• , g(x, m)f g q. 



For each t > 0, we denote by {J 7 ^, i < s < T} the natural filtration of the Brownian 
motion {W s — Wt, t < s < T}, augmented by the P-null sets of T ' . 

It follows from the results of Sections 2 and 3 that for each (t, x), there exists 
a unique triple (Y^ x , Z f ' x , K*' x ) in S 2 x M 2 x N 2 of {J*, t < s < T} progressively 
measurable processes, which solves the following RBSDE: 

Yi(s) = <?(X^(T),z) + J T 1 /,(r,X t ^(r),Y i (r),Z i (r),i)dr 



J dK t (r)- f Zi(r)dW(r), sG [0,T]; 

in| 

H 

J (Yi(s) - min{lK s ) + Hhj)}^) dK^s) = 0; i G A. 



^(s) < nm{Y-(s) + fc(iJ)}, s£ [0,T]; 



(5.1) 
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We now consider the related system of variational inequalities. Roughly speaking, 
a solution of the system of variational inequalities is a function u : [0, T] x R d — > R m 
which satisfies: 

max < —d t Ui(t,x) — £ui(t,x) — ip(t,x,Ui(t,x),Vui(t,x)a(t,x),i), 



U{(t, x) — mm(uj(t, x) + k(i, j)) > = 0, (5-2) 

j¥=i J 

(t, x, i) G [0, T] x R d x A, with the terminal condition 

Ui(T, x) = g(x, i), (x, i) eR d x A. (5.3) 
Here, the second-order partial differential operator C is given by 

More precisely, we shall consider solution of (5.2) in the viscosity sense. It will be 
convenient for the sequel to define the notion of viscosity solution in the language of 
sub- and super-jets (see, e.g., [4]). Below, S(d) will denote the set of d x d symmetric 
nonnegative matrices. 

Definition 5.1. Let u G C((0,T) x R d ;R) and (t,x) G (0,T) x R d . We denote 
by V 2,+ u(t, x) [the "parabolic superjet" of u at (t,x)] the set of triples (p,q,X) G 
Rx R d x 5(d) which are such that 

u(s,y) < u(t,x) +p(s -t) + (q,y - x) + ^(X(y -x),y- x)) + o(\s -t\ + \y- x\ 2 ). 

Similarly, we denote by / P 2, ~u(t,x) [the "parabolic subjet" of u at (t,x)J the set of 
triples (p, q, X) G R x R d x S{d) which are such that 

u(s,y) > u(t,x) +p(s -t) + (q,y - x) + ^(X(y -x),y- x)) + o(\s -t\ + \y- x\ 2 ). 

We can now give the definition of a viscosity solution of the system of variational 
inequalities (5.2) and (5.3). 

Definition 5.2. u G C([0, T] x R d ;R m ) is called a viscosity subsolution (resp., super- 
solution) of (5.2) and (5.3) ifu{(T,x) < (resp. >) g(x,i), (x,i) G R d x A, and at any 
point (t,x,i) G (0, T) x R d x A, for any (p,q,X) G V 2 ' + Ui(t, x) (resp., V 2 ~Ui(t,x)) , 

max | —p—-Tr(o~a T (t,x)X) — (b(t,x),q)—i[>(t,x,Ui(t,x),qo~(t,x),i), 

Ui(t,x) — mm(uj(t, x) + k(i,j)) > < (resp., >) 0. 
j¥=i J 

u G C([0,T] x R d ;R m ) is called a viscosity solution of (5.2) and (5.3) if it is both 
a subsolution and a supersolution of (5.2) and (5.3). 
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We now define 

Ui (t,x) :=Yl' x (t), (t, x, i) G [0, T] x R d x A; u := (m, u 2 , • ■ ■ , u m f . (5.4) 
Note that u is deterministic. 

Lemma 5.1. Xei £/ie Hypotheses 3.1, 5.1 and 5.2 hold. For each i G A, m G C([0, T] x 

Proof: For (t,x) G [0,T] x R d , i e A and o(-) G ,4j, let (C/^, V^) be the unique 
solution of the following switched BSDE: 

U{s) = g(X t ' x (T),a(T)) + A a ^(T) — A a ('\s) 

+ J ^(r,X^ x (r),U(r),V(r),a(r))dr 

-J V(r)dW(r), se[0,T]. (5.5) 

From Theorem 3.1, we have 

Ui{t,x)= inf Uf'\t), {t,x,i) G [0, T] x R d x A. 
a(-)eA\ ' 

By some classical stability arguments, we obtain the continuity of Uj. □ 

Theorem 5.1. Let the Hypotheses 3.1, 5.1 and 5.2 be true. The function u given by 
(5.4) is the viscosity solution of the system of variational inequalities (5.2) and (5.3). 

Proof: We are going to approximate RBSDE (5.1) by penalization, which was 
studied in Section 2. For each (t,x) G [0, T] x R d , let { n y*' x (s), n Z t ' x (s), t < s < T} 
denote the solution of the penalized BSDE: 

Y t (s)=g(X^(T),i) + [ T r(r,X t > x (r),Y(r),Z l (r),i)dr- f Z { (r) dW (r) , (5.6) 

J S J S 

where for (t, x, y, z, i) G [0, T] x R d x R m x R d x A, 

4) n (t,x,y,z,i) := ip(t,x,y i: z,i) - nJ^iVi ~ Vj ~ Hhj)) + - 

It is known from [15] and [1] that 

u n {t, x) := n y*' x (t), (t, x) G [0, T] x R d , 

is the viscosity solution to the following system of parabolic PDEs: 

-d t u?(t,x) -Cuf(t,x) - ^ n {t,x,u n {t,x),Wu^(t,x)a(t,x),i) = 0, 

u?(T, x) = g(x, i), (t, x, i) G [0, T] x R d x A. (5.7) 

However, from the results of Section 2, for each (t, x, i) G [0, T] x R d x A, 

u™(t, x) I Ui(t, x), as n — > 00. 
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Since u n and u are continuous, it follows from Dini's theorem that the above conver- 
gence is uniform on compacts. 

We now show that u is a subsolution of (5.2) and (5.3). Let (t,x,i) be a point in 
[0,T] x R d x A. Since u is defined by (5.4), 

Ui(t,x) < mm(uj(t,x) +k(i,j)), 

and 

Ui(T,x) = g(x,i). 

Let (t,x,i) G (0,T) x R d x A and (p,g, X) G V 2 ' + Ui(t,x). From Lemma 6.1 in [4], 
there exist sequences 

rii — > +oo, 
(fy,x/) -> (i,x), 
(Pi,9,,X,) G T^X'C*,*), 

such that 

(pi,qi,X t ) -> (p,q,X). 
On the other hand, for any Z, from (5.7), 

-Pi - 7;7>(o-cr T (*, x)X,) - (6(t, x), ft) - V(i, x, ul l (i, x),qia(t, x),i) < 0. 

Hence, taking the limit as j — > oo in the last inequality yields: 

-p ~ 7^Tr(aa T (t, x)X) - (b(t, x), q) - ip(t, x, Ui(t, x), qa(t, x),i) < 0. 

We have proved that u is a subsolution of (5.2). 

We conclude by showing that u is a supersolution of (5.2) and (5.3). Let (t,x,i) G 
(0, T) x R d x A be a point at which ui(t, x) < min^^u^i, x)+l{i, j)), and let (p, q, X) G 
V 2 '~~Ui(t,x). Again from Lemma 6.1 in [4], there exist sequences 

n\ — > +oo, 
(t h xi) -> (t,x), 
(p,,«,X,) G P 2 '-< ; (t,x), 

such that 

(pi,qi,X t ) -> (p,q,X). 
On the other hand, for any Z, from (5.7), 

- W - ^Tr(aa T (t, x)JQ) - (6(t, x), ft) - ^ n! (i, x, u ni (t, x), Ql a(t, x),i) > 0. 

From the assumption that Ui(t,x) < mmj^i(uj(t, x) + Z(i, j)) and the uniform conver- 
gence of u n , it follows that for j large enough, u™ l (ti,xi) < min,yj(u™ ! (ti, x{) + 
hence, taking the limit as j — > oo in the last inequality yields: 

-p - ^Tr(aa T (t, x)X) - (b(t, x),q) - ip(t, x, x), ftr(t, x), i) > 0. 

We have proved that u is a supersolution of (5.2) and (5.3). 

Hence, u is a viscosity solution of (5.2) and (5.3). □ 
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Remark 5.1. Uniqueness of viscosity solution to (5.2) and (5.3) follows from classical 
arguments. See, e.g. [4] and [18]. 
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